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Abstract 



We consider a binary system made of self-gravitating bodies embedded in a constant and 
uniform external field g. We analytically work out several orbital effects induced by a puta- 
tive violation of the Strong Equivalence Principle (SEP) due to g. In our calculation, wc do 
(N- not assume e ~ 0, where e is the binary's orbit eccentricity. Moreover, the spatial orienta- 

tion g of the external polarizing field is not a-priori fixed. Our results do not depend on any 
particular SEP- violating theoretical scheme. They can be applied to general astronomical 

■ and astrophysical scenarios. Finally, our findings can be used to better interpreting present 
I and future SEP experiments, especially when several theoretical SEP-inducing mechanisms 

■ may be at work, and to suitably designing new dedicated tests as well. 

o: 

CN| ■ Keywords: Experimental studies of gravity; Experimental tests of gravitational theories 

PACS: 04.80.-y; 04.80.Cc 

K*^ ! 1 Introduction 

H 



If bodies with non-negligible gravitational binding self-energies, like astronomical and astro- 
physical objects (planets and their natural satellites, main-sequence stars, white dwarfs, neutron 
stars), move with different accelerations in a given external field, the so-called Strong Equiv- 
alence Principle (SEP) would be violated. While the General Theory of Relativity (GTR) 
assumes the validity of SEP, it is generally violated in alternative theories of gravity; see [H [2] 
and references therein. 

Potential SEP violations in the Earth-Moon system freely falling in the external gravitational 
field of the Sun, theoretically predicted by Nordtvedt long ago OS], are currently searched for 
with the Lunar Laser Ranging (LLR) technique [5l [H [7] in the framework of the Parametrized 
Post-Newtonian (PPN) formalism. There are projects to test SEP with the Martian moon 
Phobos and the Planetary Laser Ranging (PLR) technique j8] as well. 

As shown by Damour and Schaefer [9] , the acceleration due to SEP violation occurring for a 
two-body system in an external gravitational field g like, e.g., the Galactic field at the location 
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of a binary pulsar, is, to leading ordei0, 

AsEP = (Ap- Ac)£/, (1) 

where p and c denote the pulsar and its less compact companion; for a generic body with non- 
negligible gravitational self-energ}{l f g^av j A accounts for the SEP- violating difference between 
the inertial mass and the gravitational mass m^: at first post-Newtonian level, it is 

A — — ^ — n c c — ^g''^^ (o\ 

mi miC^ 

where c is the speed of light in vacuum. The external field g in eq. ([T|) is assumed here to 
be constant and uniform, so that it induces a gravitational analog of the Stark effect [9]. To 
maximize the size of potential SEP violations, one should look at compact objects like neutron 
stars: indeed, it is 



^grav 



1.9 X 10-", |e®av| ~ 4.2 X IQ-^o, \e%,,\ ~ 1.3 x 10"^, \ef^^\ ~ 10"!. (3) 



The SEP-violating parameter r]i can be parameterized in various ways depending on the the- 
oretical framework adopted; see [H [2] for recent overviews. As far as the Earth-Moon system 
and the PPN framework is concerned, latest published bounds on the Nordtvedt parameter are 

7?i = ?7N = (4.0±4.3) X 10"^ [6], (4) 



= r?N = (-0.6±5.2) X 10"^ [7|. (5) 

Williams et al. in a recent analysis [13| released 

r/i =r/N = (-1.8±2.9) X 10"^. (6) 

In this paper, we deal with the orbital effects due to eq. ([T|) in a uniform way by obtaining 
explicit and transparent analytical expressions valid for the main observable quantities routinely 
used in empirical studies. In Section [2] we work out the SEP rates of change of the osculating 
Keplerian orbital elements, which are commonly used in solar system and binaries studies. 
The SEP-like shift of the projection of the binary's orbit onto the line-of-sight, which is the 
basic observable in pulsar timing, is treated in Section [3l In Section S] we calculate the SEP 
perturbation of the radial velocity, which is one of the standard observable in spectroscopic 
studies of binaries: in systems suitable for testing SEP, the pulsar's companion is a white dwarf 
or a main sequence star for which radial velocity curves may be obtained as well. The primary- 
to-companion range and range-rate SEP perturbations, potentially occurring in systems like 
the Earth-Moon one, are calculated in Section [5l Section [6] is devoted to summarizing and 
discussing our results. We will not resort to a-priori simplifying assumptions about the binary's 
orbital eccentricity. Moreover, we will not restrict ourselves to any specific spatial orientation 
for g. Our results are, thus, quite general. They can be used for better understanding and 
interpreting present and future SEP experiments, hopefully helping in designing new tests as 

^It is assumed that the external field g is weak at the location of the binary system considered, and that its 
sources are mainly made of non strongly self-gravitating bodies [2]. 

■^For a non-relativistic spherical body of uniform density, it can be cast [lOIIll) Sgrav ~ 3M^G/5R, where G is 
the Newtonian constant of gravitation, M is the body's mass, and R is its radius. For highly relativistic neutron 
stars, see, e.g., [12]. 
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well. Finally, we remark that our calculations are model-independent in the sense that they 
do not depend on any specific theoretical mechanism yielding SEP violations. Thus, they may 
be useful when different concurring SEP-inducing scenarios are considered like, e.g., MOND 
and its External Field Effect [14^ [T5] , Galileon-based theories with the Vainshtein mechanism 
|16[ I17j. and variations of fundamental coupling constants as well |18 [ I19 t [2U]. Moreover, there 
are also other effects, like Lorentz symmetry violations [21], which affect the same SEP peculiar 
observables like the eccentricity |22| . 



2 The long-term rates of change of the osculating Keplerian 
orbital elements 

The SEP acceleration of eq. ([1]) is much smaller that the standard two-body Newtonian and 
Post-Newtonian ones; thus, its orbital effects can be treated perturbatively. It can be thought 
as due to the following perturbing potential 

U = E g-r, (7) 

where 

H = Ae - Ap : (8) 

r is the relative position vector of the binary. 

The Lagrange planetary equations [23] allow to work out all the orbital effects caused by 
eq. d?]) in an effective way which encompasses just one integration. 

Damour and Schaefer [9] used a different formalism and an approach requiring six inde- 
pendentH integrations; moreover, they did not release complete explicit expressions written in 
terms of the standard Keplerian orbital elements, and only the rate of change of the eccentricity 
was considered in detail. Freire et al. [2], using the same formalism as Damour and Schaefer 
|S], dealt with the inclination as well. 

As a first step, the perturbing potential of eq. ([7]) must be evaluated onto a reference 
trajectory assumed as unperturbed with respect to the effect we are interested in. To this aim, 
we will adopt the standard Keplerian ellipse, although it would be possible, in principle, to 
use a Post-Newtonian orbit as well [Ml [25] • We do not use such a relativistic unperturbed 
path since it would only yield additional SEP-IPN cross terms which are quite smaller than 
the already small SEP-OPN effects. Then, an average over one full orbital period of the binary 
must be taken; by using the eccentric anomaly E as fast variable of integration, the result is 

(U) = ^— ^ {cos a; {gx cos Q + Qy sin Q) + sinuj [g^ sin / -|- cos / {gy cos Q — cjx sin J?)]} , (9) 

where g = \g\ ,g = g/g, and a, e, /, w, Q are the semimajor axis, the eccentricity, the inclination, 
the argument of pericenter and the longitude of the ascending node of the binary orbit referred 
to a generic coordinate system. As far as binary pulsar systems are concerned, the reference 
{X, y} plane is the plane of the sky tangent to the celestial sphere at the location of the system 
considered. The inclination I of the binary's orbital plane refers just to the plane of the sky, 
so that the reference Z axis is directed from the solar system's barycenter to the binary along 
the line-of-sight direction; the reference Y axis is usually directed towards the North Celestial 



■^Averages of the time derivatives of the orbital energy, the orbital angular momentum and the Laplace- 
Runge-Lenz vector were taken in |9]- The orbital angular momentum and the Laplace- Runge-Lenz vector are 
perpendicular. 
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Pole, and the reference X axis, from which the node H is counted, is in the east direction. The 
unit vectors Io,Jo,Kq of such a coordinate system can be expressed in terms of the Celestial 
right ascension (RA) a and declination (DEC) 6 of the binary as |26[ f27\ [25] 



— sma, 
cos Q, 
0, 



(10) 



— sin 5 cos a, 

— sin (5 sin a, 
cos 5, 

cos 5 cos a, 
cos 5 sin a, 
sin 6. 



ill) 



(12) 



Instead, in studies pertaining our solar system the reference {X, Y} plane is customarily chosen 
as coincident with the mean equator at J2000.0, so that the reference X axis is directed towards 
the Vernal Equinox Y at J2000.0 and the Z axis points towards the North Celestial Pole. In 
obtaining eq. Q we have reasonably assumed that g does not vary over a full orbital revolution; 
moreover, we did not make any a priori simplifying assumption about the spatial orientation of 
g which, in general, will not be directed along any particular direction of the coordinate system 
adopted. In the case of a binary pulsar system, the SEP polarizing field g can be assumed to 
be approximately directed towards the Galactic Center (GC), whose Celestial coordinates are 
[29] 



ace = 17''45°'37'.224 = 266.4051 deg. 



(13) 



5qc = -28°56 10 .23 = -28.936175 deg. 



(14) 



Thus, by projecting g onto eq. ()10p -eq. (]12p . it can be expressed in terms of a, 6 of the binary 
and of OGC^GC as 



cos Sgc sin {a - ace) , 

cos (a — ogc ) cos 5gc sin 5 + cos 5 sin Jqc , 

I + cos (a — ace) cos (5gc cos 5 + sin 5 sin 5gc- 



(15) 



The Lagrange planetary equations [23] and eq. (Hj) straightforwardly yield the following 
long-term rates of changes of the six standard osculating Keplerian orbital element^, and of 



*M = rib (t — to) is the mean anomaly, where to is the time of passage at pericenter and rib = a/ GM/a^ is 
the Keplerian mean motion. 
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the longitude of the pericenteiH vj as well. 

da , 



(16) 



3H<7^/r 



[cjz sin / cos uj + cos / cos oj {gy cos Q — Qx sin Q) 



— sintj {qx cos n + cjy sin /?)] , 

dl \ SHf^e cos w [qz cos / + sin / {cjx sin fl — Qy cos 
/ ~ 2nbaVl - 



(17) 
(18) 



dQ \ 3H(j(e CSC I sin u; [^^ cos / + sin / {cjx sin Q — gy cos J?)] 
dt / 2nba\/l - 



(19) 



(ii / 2enba\/l - 



{ (— 1 + e^) cos oj {cjx cos Q + (jy sin i?) + 



+ sin [— 5y cos / cos /2 + (e^ esc / — sin /) + gx cos / sin Q\ } , 



(20) 



dvj 



dt I 2en\,a\/\ — & 



{ (l — e^) cos OJ (ijx cos H + tjy sin f?) + 



+ sinw 



(cos / — e^) (g^ cos Q — Cjx sin j?) + ^2(1 — 6^ + cos /) t 



an 



(21) 



dM 
"df 



3Hg (1 + 
2enba 



{cos u {gx cos Q + (jy sin J7) 



+ sin UJ [gz sin / + cos / {(jy cos i7 — gx sin /?)]}. 



(22) 



About the rate of the mean anomaly, we recall that in obtaining it perturbatively the mean 
motion rib has to be considered constant; thus, eq. ([22]) is proportional to the rate of the time 
of pericenter passage. In obtaining eq. (|16|) -eq. (|22|) . we did not make any a-priori simplifying 
assumption on the orbital configuration, i.e. we did not assume e — )• 0, / — )• 0. In the limit 
e — )• the orbital plane remains unchanged since the long-term changes of both the inclination 
/ and the node Q vanish. Instead, the rate of change of the eccentricity e remains substantially 
unaffected since it is not of order 0{e). In the limit of small eccentricities and small inclinations, 
the mean longitude I = w + Ai is customarily used. From eq. (|2ip -eq. (j22p it turns out that 



^It is a "dogleg" angle defined as vj = Q -\- oj. It is used in some specific cases such as, e.g., analyses of 
planetary motions in our solar system. 
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its rate is of order 0{e): 



' dl \ 3"6 

— ) — rib = — {3 cos a; {gx cos Q + Qy sin Q) + sin a; [cjy (2 + cos/) cos Q- 

dt / 4nbi3 



'/ 

Qx (2 + cos /) sin Q + (jz sin / + 2gz tan ( — 



(23) 



Thus, in the small eccentricity limit, all the osculating Keplerian orbital elements remains 
unaffected by a SEP violation, apart from the eccentricity. Note that eq. (jl6p -eq. (j22p and 
eq. ([23]) are all proportional to n^^a~^ = a^/^, so that they are larger for detached binaries. 
For large eccentricities the node and the inclination variations can, in principle, be used as well 
as further observables. Interestingly, eq. (jlSh -eq. ()19p show that their ratio is independent of 
both g and g] moreover, it depends neither on a nor on e, nor on Q amounting to 



esc/ tan tj. (24) 



/ 

For systems exhibiting large inclinations, eq. (j24p may have relevant observational implications 
thanks to its peculiar signature, provided that both the inclination and the node are accessible 
to observations. Also the ratio of eq. (jlSp and eq. ()17p could be considered [2]; however, if on 
the one hand the ratio of the rates of the inclination and the eccentricity does not depend on g 
and a, on the other hand it depends on g and on e, /, i7, uj as well, contrary to eq. (j24p which 
is simpler. 

Damour and Schaefer [S] define the following vector directed along the external fielcjfl 

Then, they use the unit vector a directed along the semimajor axis towards the pericenter, the 
unit vector k directed along the orbital angular momentum, and b = k X a. Their expression 
for the averaged rate of change of the eccentricity is [9] 

de \ 



9) = VT3^(6./). (20) 



A direct comparison of eq. ()26l) with our eq. ()17p can actually be done by expressing the unit 
vectors k and d in terms of the Keplerian orbital elements as 



k 



sin /sin 12, 

— sin /cos /?, (27) 
cos /. 



and 



cos n cos u — COS / sin il sin w, 
d = sin i? cosw + COS /cos /? sincj, (28) 
sin/ sin tj. 



^Note that the quantity / has the physical dimensions of reciprocal time, i.e. [/] — T 



6 



By using eq. ([27j) -eq. ([28]) . it turns out that eq. ([26]) coincides with our eq. pT|) . 

As far as the incHnation I is concerned, Freire et al. [2], by adopting the same formahsm of 
Damour and Schaefer [S], obtain 



(j28p it can be shown that eq. ()29p agrees with our eq. (jlSp . 

Neither Damour and Schaefer nor Freire et al. [2] exphcitly considered the SEP-induced 
rates of change of the node (eq. p^ ). the pericenter (eq. ([20]) -eq. (pT]) ) and the time of 
pericenter passage (eq. ([2^ 1. 

In binaries hosting a radiopulsar, the rate of change of the inchnation can be expressed in 
terms of the rate of change of the pulsar's projected semimajor axis Xp = Op sin//c, where c is 
the speed of light in vacuum, as / = (ip/xp) tan /. About the measurability of the node Q in 
wide binaries hosting a radiopulsar, it could be determined only if they are close enough to the 
Earth. Indeed, in this case the orbital motion of the Earth changes the apparent inclination 
angle / of the pulsar orbit on the sky, an effect known as the annual-orbital parallax [26]. It 
causes a periodic change of the projected semi-major axis. There is also a second contribution 
due to the transverse motion in the plane of the sky [27], yielding a secular variation of the 
projected semi-major axis. By including both these effects in the model of the pulse arrival 
times, Q can be determined, as in the case of PSR J0437-4715 [28], located at only 140 pc 
from us. Interestingly, such an approach may also be used in optical spectral observations 
of binary stars possessing a sufficiently well determined radial velocity curve [27]. As far as 
exoplanets are concerned, spectroscopic variations of the hosting star during the transits can, 
in principle, be used to measure Q [30]. According to Fluri and Berdyugina [31], the analysis 
of the polarization of the light scattered by the planetary atmospheres may allow to determine, 
among other things, fl as well. 

The SEP orbital rates of eq. (116p -eq. (I23p can be defined as long-term effects since they were 
obtained by taking averages over one full orbital period Pb of the binary system. It must be 
noted that, strictly speaking, they cannot be considered as secular rates. Indeed, in calculating 
them, it was assumed that the external field g was constant over Pb- The extent to which 
such an approximation can be considered valid depends on the specific system considered. It 
is certainly true for binary pulsar systems, for which timing data span some decades at most, 
when g is due to the Galaxy: in this case, g can reasonably be considered as constant over any 
foreseeable data analysis based on existing (and future) pulsar timing records because of the 
extremely low variation of the Galactic field at the pulsar's location during such time spans. 
If, instead, g is due to a remote third body with a small orbital frequency n^^ with respect to 
TT-b, a relatively low modulation may be introduced by a slowly varying g{t). A quite different 
situation occurs for systems like the Earth-Moon one whose external field g is due to a relatively 
close body like the Sun having an orbital frequency n^^ comparable to rib. In this case, g cannot 
be considered as constant, and a second integration has to be performed over the orbital period 

of the third body. 

In eq. ()16p -eq. (j23p the Keplerian orbital elements /, /?, uj determining the orientation of the 
binary orbit in space are present. They are fixed only in the two-body Keplerian problem for two 
pointlike masses. In any realistic situation, they are slowly varying with frequencies typically 
far smaller than nb. It is true even if no other bodies are part of the system because of general 
relativity itself causing the well known IPN gravitoelectric and gravitomagnetic precessions. 




(29) 



where Kq is the unit vector of the line-of-sight. Also in this case, with the aid of eq. (|27p -eq. 
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Thus, even if g can be considered constant, in general it is not so for /, n,uj. Again, the 
pecuharities of the specific binaries at hand may greatly reduce the variability of them. For 
example, the more the system's orbit is large, the weaker are the effects of general relativity. 

Our results of eq. (fT6|) -eq. ([23|) are quite general, and can be applied to any gravitationally 
bound two-body system immersed in a (weak) external field: they are valid for any orbital ge- 
ometry of the binary system and for a quite generic spatial orientation of in a given coordinate 
system. 



3 The line-of-sight perturbation Ap 

In timing of binary systems typically hosting radiopulsars, the basic observable quantity is the 
projection p of r along the line-of-sight [S] since its variations Ap are straightforwardly related 
to the timing via At = Ap/c. In view of the fact that p = Z, one has to consider the Z 
component of the binary position vector r 



Z = r sin / sin (iv + i/) , 



(30) 



where ly is the true anomaly. The SEP perturbation of p due to eq. ([T]) can be straightforwardly 
obtained in the following way: 



(ApsEp) 
where 



fdZ\ , 



2tt 



dZ dE dM ^p dZ dn 



dEdM dt 



a,e,I,uj, (31) 



a (1 — ecosE) 



(32) 



cos z/ 



COS E — e 
1 — e COS E ' 



(33) 



sm u 



\/l — sinE 
1 — e cos E ' 



(34) 



1 



dE 



dM l-ecos£:' 



(35) 



dt 1 — e cos E 
dE Tib ' 



(36) 



and dn/dt are the instantaneous rates of change of the Keplerian orbital elements computed 
onto the unperturbed Keplerian orbital ellipse according to, e.g., the right-hand-sides of the 
Gauss perturbative equation^ [23] . 



^In principle, dM/dt includes also Uh, which yields the purely Keplerian part (ApKop) of the shift of the 
line-of-sight component of the orbit. It vanishes. 
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We have 



( A/)sEp) = ^ ^ A 2 { ~ + 3'^^) sin / sin 2w+ 

+ (4 — 6e^ + 3e^) cos / sin 2w (5^ cos Q — Qx sin /?) + (37) 

+ [— + (4 — 6e^ + 3e'*) cos 2w] (^a; cos il + gy sin /?) | . 

No approximations in and / were used in deriving eq. (|37|) . A purely formal singularity 
appears in eq. ([37]) for e — )• 0: actually, it is unphysical since it is cured by using properly 
chosen non-singular orbital elements (see, e.g., [32)). 

4 The radial velocity perturbation Ap 

A standard observable in spectroscopic studies of binaries is the radial velocity. Up to the 
velocity of the binary's center of mass Vq, it is, by definition, dp/dt, i.e. the temporal rate of 
change of the line-of-sight projection of r. 

The unperturbed, Keplerian expression for p can straightforwardly be obtained from eq. 
([30]) and eq. (gOD as 

dZ du n\,a sin/ 

^ = 'd^dM'''' = [ecosu; + cos + z.)] . (38) 

The SEP-induced perturbation (psEp) on p can be worked out by replacing Z with eq. ()38p in 
eq. ([3T]) . For computational purposes, it turns out more convenient to use the true anomaly u 
as fast variable of integration instead of the eccentric anomaly E. 
By using 



1 + e cos ly ' 
di^ / a \ 2 



(39) 



dM \r 
dt {l-e^f^ 



Vl-e2, (40) 



di^ nb (1 + ecos v)^ ' 

one finally gets 



(41) 



/A ■ \ 97r(7H\/l — sin/ ,^ . ^ 

( ApsEp) = 5 {9z sm / cos 2uj+ 

e^nb 

+ cos / cos 2oj {cjy cos Q — cjx sin /?) — sin 2a; [g^ cos Q + gy sin J7)} . (42) 

Also in this case, the use of suitably defined non-singular orbital elements cures the formal 
singularity in eq. (I42p for e — t- 0. 
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5 The range and range-rate perturbations Ar and Af 



By proceeding as in the previous Sections, it is possible to work out the SEP-induced range 
and range-rate perturbations Ar and Ar, respectively. More precisely, we used eq. ([31]) in 
which we replaced Z with eq. (|39p for Ar, and with dr/dt = n\^aesm.v / y/\ — for Ar. Such 
observable quantities are routinely measured in, e.g., the Earth-Moon LLR experiment and in 
several spacecraft-based interplanetary missions. The following calculations refer to the vector 
r connecting the centers of mass of the two bodies constituting the binary system; thus, they are 
valid when the distances of the ranging devices (laser retroreflectors, transponders on probes, 
etc.) from the centers of mass of the bodies are negligible with respect to r. 
We obtain 

SttH-v/I — 

(ArsEp) = 2 ^ ^™ ^ ^ '^^^ ^ ~ 9x sin /?) — 

— sincj {qx cos Q + cjy sin /?)] , (43) 

,^ . , 37rH[l + e(2-e)] , ^ ^ ■ r^^ 

(ArsEp) = . _ 2\ ^'^^^ ^ \3x cos il + Qy sm 12) + 

-|- sin a; [g^ sin / -|- cos / (gy cos Q — g^ sin /?)]} . (44) 
No approximations in e were used in deriving eq. (|43p and eq. (|44p . 



6 Summary and conclusions 

We systematically worked out all the relevant orbital effects induced by a Stark-like SEP viola- 
tion in a binary system made of two self-gravitating bodies immersed in an external polarizing 
field assumed constant and uniform with respect to the characteristic temporal and spatial 
scales of the binary. 

We provided the reader with full analytical expressions for the SEP-induced rates of change 
of all the six osculating Keplerian orbital elements (Section [2]), for the projection of the binary's 
orbit along the line-of-sight (Section [3]) and its time derivative (Section U]), and for the range 
and range-rate (Section [5|). We did not make any a-priori assumption about the orientation of 
g in space. We did not make simplifying assumptions about the orbital geometry of the binary. 
Thus, our calculations are quite general, and are valid, in principle, for all astronomical and 
astrophysical scenarios. 

It is important to have at disposal explicit expressions of the different relevant SEP effects 
since each one has a peculiar temporal pattern which may be helpful in separating it from other 
possible competing signatures. For example, PSR B1620-26 is a pulsar-white dwarf binary 
in a relatively wide orbit (Pb ~ 2 x 10^ d) orbited by quite distant circumbinary planet-like 
companiorEI (/^ ~ 3 x 10^ d) [36l EZ]. It is not possible to treat here the impact of other 
dynamical effects acting as unwanted distubances as far as SEP violations are concerned since, 
in general, they vary from a system to another. Such a kind of analysis can be done with the 
specific scenario at hand. For a discussion of the rate of change of the eccentricity in binaries 
hosting compact objects, see, e.g., [2]. 

^For exoplanets, see [33l [Ml [35] . 
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Our results are important also to accurately asses the overall uncertainty which can be 
obtained when constraints on A are inferred by comparing them with the corresponding em- 
pirically determined quantities. Indeed, taking into account only the accuracy with which an 
observable like, say, de/dt can be determined is, in principle, not enough; a propagation of 
the errors affecting all the parameters, like g = V"Jj,p/d [351 ES]) accurate to roughly ~ 10% 
according to the uncertainties in the Galactic circular velocity Vdr [lO], and its orientation g, 
entering the corresponding theoretical prediction must be done as well to correctly infer the 
total, systematic uncertainty in A. 

Our expressions are also useful in designing suitable SEP tests and in better interpreting 
present and future experiments, especially when different competing theoretical mechanisms 
yielding SEP violations are considered. Indeed, we found that the ratio of the node and incli- 
nation SEP-type precessions depends only on the inclination itself and on the pericenter, which 
is a quite distinctive signature. 
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